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REPORT 



DYNAMIC LOADS ON AIRPLANE STRUCTURES DURING LANDING 
By M. A. Blot and R. L. Bisplinghoff 



1. SUMMARY 



The application of transient theory to the rational 
determination of dynamic loads on airplane structures during 
initial landing impact i8 discussed. Simplified procedures 
are described by which the distribution of the maximum at- 
tainable value of the dynamic stresses in the structure may 
be obtained. Illustrations of the procedure are given by 
numerical examples for the case of airplane wings. This 
indicates approximate orders of magnitude to be expected in 
a typical problem. The validity of the underlying simplifi- 
cations and assumptions is discussed. A brief outline of 
the general mathematical theory o-f transients in undamped 
elastic system is presented. 



2. INTRODUCTION 



During landing, the airplane structure undergoes 
transient oscillations which are excited by the initial land- 
ing impact , Recent experience has shown that in the case of 
large aircraft these oscillations may produce critical design 
conditions for the structure; whereas heretofore design loads 
for the landing condition have been based upon calculations 
which assume the structure to be rigid. Since the advent of 
larger aircraft has resulted in heavier and more flexible 
structures, these calculations are considerably in error. 
This together with the fact that flight load factors are re- 
duced in the case of large aircraft makes it necessary that 
methods be developed for predicting dynamic landing loads. 
The present theoretical investigation of this problem is of a 
preliminary nature. It has the twofold purpose first of pre- 
senting methods by which the designer may predict the dynamic 
loads and second of serving as a guide in the experimental 
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investigation "by defining the significant factors involved 
and determining the approximate magnitude £f the quantities 
to be measured. Wk'T': '* • " * 5 

It appears that the problem may be approached in two 
different ways: 

In the case of the landplane to consider the airplane 
structure and its landing gear as a whole, and to introduce 
the actual force displacement characteristics of the landing 
into the theory. In this procedure the dynamic stresses re- 
sult from the sudden application of moving constraints im- 
posed on the airplane during landing. Similarly for the 
seaplane the elastic structure and the water surrounding the 
hull may be considered as interacting bodies. While this 
method is not precluded in the investigation of specific 
cases or for research purposes, it involves inherent complex- 
ities, such as those resulting from the nonlinear properties 
of the landing gear and the variable-mass effects of the 
water surrounding the hull of a seaplane, which tends to make 
t.his type of approach less adequate for design purpose. 

In the other procedure, in which stresses in the struc- 
ture are considered to be caused by a landing impact force 
appliod directly to the structure, it is assumed that the 
time history of the impact force may be investigated inde- 
pendently of the elastic properties of the structure. In 
this way the investigation involves two separate phases - a 
study of the landing forces and a study of the dynamic be- 
havior of the structure under such forces. This procedure 
involves the assumption that a landing impact force may be 
defined in such a way that its time history is for all prac- 
tical purposes independent of the elastic properties of the 
structure. Since it is believed that the relative simplicity 
of the latter approach overshadows the approximations in- 
volved, it has been adopted as the basis for the present work. 
This procedure was described previously in reference 1 in 
connection with the problem of dynamic stresses in buildings 
during an earthquake and the present work is essentially an 
adaptation to aircraft structures of the methods described 
in this reference. 

It is assumed that in first approximation the damping 
and the aerodynamic forces may be neglected. The landing 
impact force is applied directly to the elastic structure as 
an external force of given time history. The response of 
the structure is represented as a superposition of natural 
modes excited by the landing impact. A first step in the 
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analysis is therefore to obtain the natural modes of the air- 
plane either by calculation or by a shake test. Calculation 
methods have recently been developed by which natural modes 
of airplanes may be derived in a relatively simple way. 

An important feature of the designer's approach to the 
landing loads problem is the fact that he is not so much in- 
terested in the actual time history of the structure as he 
is in the magnitude of the highest attainable stresses during 
the operation of the airplane. This viewpoint was introduced 
in the procedtire by using a statistical approach. The stress 
amplitudes of each mode are superposed with their positive or 
negative values irrespective of phase and the worst possible 
combination is used as a basis for design. Furthermore, the 
stress history in each mode is not actually computed but the 
stress amplitude is obtained directly from a graph represent- 
ing what is designated as a "dynamic response factor." This 
factor itself results from a statistical analysis of the ef- 
fect of forces of various time histories on a single degree 
of freedom oscillator, vising a sufficient number of such time 
histories to represent all possible types of landing condi- 
tions. Values of the dynamic response factors are obtained 
by applying typical time history excitations to a torsional 
pendulum (described in reference l) and measuring the maxi- 
mum amplitude of its response. In this way, when the natural 
modes of the airplane are known, the landing loads are read- 
ily obtained without the necessity of integrating the dif- 
ferential equations which govern the behavior of the elastic 
structure in the transient condition. 

The method has its limitations and is not necessarily 
applicable to all types of transient problems. Some of 
these limitations are pointed out in the di ecus si on, but the 
extent to which the methods are valid for some particular 
class of problems will have to be determined experimentally. 



3. NOTATION 



M,m mass 

k spring constant 

Q,(t) generalized force 

w circular frequency 
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T period 

t,T time variable 

V dynamic response factor 

q generalized coordinate 

z deformation displacement of any point on the wing 

$ normal function describing the wing mode shape 

a normal function describing the mode shape of twisting 
about the clastic axis 

h normal function describing the mode shape of bending 

of the olastic.axis 

I moment of inertia 

S static mass moment 

j subscript denoting jth mode 

k subscript denoting kth spanwise wing station 

p mass per unit volume of prismatic beam 

A cross- sect i onal area of prismatic beam 

E modulus of elasticity 

n vertical landing load factor 

3f(t) landing reaction 

N(t) shock strut axial force - time relation 

D(t) drag force - time relation caused by wheel spin-up 

a(t) observed acceleration time history in a drop or 
landing test 

W gross weight of airplane 

W, gross weight of airplane less landing gear weight 
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effective rolling radius of wheel 




landing speed 



coefficient of sliding friction of tire on runway 
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str ess 



I 



moment of inertia of wheel and tire about axis of 
rotation 




period of shock strut axial impulse or seaplane 
vertical impulse 




period of drag impulse 
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angular displacement of wheel 



4. THE EVALUATION OF TRANSIENT MOTION OF ELASTIC BODIES 



As pointed out in the introduction, the theory proceeds 
on the assumption that the landing impact force is known. 
In this way the theoretical problem is reduced to the evalu- 
ation of the response of an elastic structure to a force of 
given time history. Methods for the determination of this ■ 
landing impact force from test data will be discussed later 
in section 5. 

In order to introduce the fundamental concepts involved 
in the present theory, the simplest possible elastic system 
will first be considered. This system is illustrated in fig- 
ure 1 as a single mass oscillator. 



3Y THE USE OF GENERALIZED COORDINATES 




A 




Figure 1 
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Denoting "by M the mass and u; the natural frequency of 
oscillation (oj is the circular frequency expressed in 
radians/sec), the spring constant Is k = M o> a . Tlie expres- 
sion giving the displacement q of this mass under a force 
F(t) of arbitrary time history is well known (references 2, 
3, and 4). It may be written as: 



where T is a variable of integration. It is usually desig- 
nated as Duhamel's integral. According to this formula the 
computation of the displacement q at the instant t re- 
quires the evaluation of a definite integral between the 
limits of integration 0 and t, and the time history of q 
is obtained by repeating this process for every value of t. 
It is noted that even in the simple case of a single-mass 
system the process of computing the transient response is 
quite elaborate. Fortunately, this difficulty may be avoided 
in adapting the theory to practical problems of design by 
considering only the maximum value of q instead of its com- 
plete time history. How this is achieved will be explained 
later (sec. 7). 

Assume for the present that the complete time history 
of the deformation in the airplane structure is desired. 
Such a structure differs from the simple system of figure 1 
by two features: 

(a) It is a free body. 

(b) It is an elastic body with many degrees of freedom. 

In order to show how the previous discussion may be ex- 
tended to include these features, consider a free system of 
two masses M and m connected by a spring of constant k 
as illustrated by figure 2. 



t 




T )d t 



(1) 



m 




M -4 — r(t) 
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Figure 2 
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The mass M is under the action of a force F(t) of arbi- 
trary time history. The straightforward way of describing 
the motion of this system is in terms of the displacements 
Xi and x 3 of each mass. However, there is a more general 
approach. By an elementary analysis it is shown in appendix 
II that this motion may be described as the superposition of 
two configurations. One in which the two masses move to- 
gether as a rigid body (fig. 3a), the other in which the cen- 
ter of gravity remains fixed with the masses moving in oppo- 
site phase and with amplitudes inversely proportional to the 
respective masses (fig. 3b). 



<lo 



m 



It 



3i 



s / / / / / s ; — / 7 f '"■> r 
(a) Higid configuration 



Figure 3 



Fixed e.g. 




7~7 / / / 7 — / / /■ 
(b) Motion about e.g. 



7^ 



Each of these configurations has only one degree of freedom. 
The displacement in the first configuration is measured by 
the quantity q 0 and in the other by the quantity q x . The 

actual . di splaceroent s of the masses M and m in terms of 
the motion of each configuration are, respectively, 

*i - <lo + I *i 

(2) 

x 2 = q 0 ~ <Li 



These two configurations may be interpreted physically as 
representing the natural modes of oscillations of the system. 
Prom this viewpoint the rigid translation q 0 may be thought 

of as a mode of zero frequency or "zero mode." The other 
configuration where the center of gravity remains fixed and 
the masses move in opposite phase represents a mode of fre- 
quency u>i. Since the coordinate q 0 of the zero mode rep- 
resents the motion of the center of gravity, it is determined 
by the motion of a single mass M 0 = M + m under the force 
E(t) (fig. 4a). As shown in appendix II, the motion in' the 
mode defined by ^ may be determined from that of an equiv- 
alent system which is constituted of a single mass M a 
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elastically restrained and under the action of a force Qj.(t) 
proportional to P(t) (fig. 4b). 



-r(t) 



(a)' 



Pilars 4 




©l(t) 



The natural frequency of this equivalent system is the same 
COi as the natural frequency of the mode which it represents 
in the actual physical system. Such coordinates as q 0 and 

q x are called generalized coordinates . The mass Mj, of the 
equivalent system is referred to as the generalized mass of 
the corresponding mode and Q x (t) is referred to as the 
generalized force for this same mode. It is also shown (see 
appendix II) that the value of this generalized mass is de- 
rived quite simply by expressing that the kinetic energy T 
in the equivalent system is the same as in the corresponding 
mode 



T = I M x q\ s = I m qi * + I M <J> 3 qi 3 
2 2 2 



(3) 



where <t> = — is the ratio of the amplitudes of the masses m 
M 

and M in the q, mode. 



Hence 



M x = m + M <t> c 



(4) 



Similarly the generalized force Q x (t) is derived by 

expressing that the work done "by F(t) on the actual system 
is equal to the work done by fti(t') on the equivalent sys- 
tem - that i s , 



or 



<t> q x P(t) = Q,x(t) <i x 



^(t) = <t> y(t) 



(5) 



The problem of finding the transient motion of two masses 

connected elastically has thus been reduced to that of two 

independent single masses for which the motion may be ex- 
pressed by Puhamel ' s integral (equation (l)). 
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Consider now a simplified airplane in which the wing is 
represented "by two masses elaetically connected to a rigid 
fuselage as illustrated in figure 5a. 




Obviously, for symmetric deformations this simplified 
structure is equivalent to the two-mass free system of fig- 
ure 2. The motion under the landing impact force may "be 
described as the superposition of a rigid translation (fig. 
5b) which represents the motion of the center of gravity and 
a natural mode (fig. 5c) which represents the motion about 
the center of gravity. The equivalent single mass systems 
are the same as in figure 4. 

This procedure may be immediately generalized to a com- 
plex airplane structure. Instead of a single deformation 
mode as in the case of the simplified airplane discussed 
above, there are actually an infinite number. It can be 
shown that the deformation of the structure may be repre- 
sented by a superposition of these modes. If damping is 
neglected as is the case here, these modes are uncoupled. 
The behavior of each one under the landing impact force is 
independent of the other. The motion of each mode is repre- 
sented by that of an equivalent single mass oscillator of 
the same natural frequency and excited by a generalized 
force. As an example, take the case of a wing in pure bend- 
ing. There are an infinite number of bending modes, a few 
of which are represented in figure 6. The motion of the 
center of gravity represented by the rigid airplane with the 
generalized coordinate q 0 . 18 referred to in the present 
text as the "zero mode." 
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'Figure 6 
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The amplitudes of the wing tip are usually selected as 
generalized coordinate and denoted "by qi for the first 
mode, q s for the second mode, and so forth. The correspond- 
ing generalized masses are determined "by the condition that 
the kinetic energy he the same for the mode and its equivalent 
system. This yields the expression: 



k 



M 3 = ^jj^k^ 2m k» 80 forth 



(6) 



In these expressions, ^ ^ 



(J) 



represents tho amplitude at 
station k of the mass m k of that station in the jth mode 
for a unit deflection of the wing tip. Similarly, the gen- 
eralized forces determined by the conditions that the work 
done by F(t) he the same in the particular mode as the 
work done hy Q(t) in the equivalent system are 



0 (t) = F(t) <t>( 3 \ and so forth 

'3 J? 



(7) 



The subscript F refers to the station at which the landing 
impulse force is applied. 
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5. THE LANDING IMPACT FORCE F(t) 



In the preceding discussion, the time history F(t) of 
the externally applied landing impact force is assumed to "be 
known. In the case of a landplane, the forces are applied 
to the wing through an oleo 6trut. In the case of a flying 
"boat, they are transmitted to the wing through a hull or 
pontoon. 

For the landplane, a convenient source of information 
of axial strut characteristics is the oleo drop test. In 
the drop test, a mass-oleo strut system is dropped in a Jig, 
and the acceleration time history of the mass is measured. 
The mass corresponds to the zero mode mass. 



An additional source of landing force data is from 
actual landings with acceler ometers installed in the airplane 
which are capable of recording time history. This method 
provides the only present source of seaplane data. A sketch 
of such a record is 6hown in figure 7. 



The high-frequency oscillations must he disregarded 
since they represent oscillations of the structure relative 
to its center of gravity. The average dotted line which is 
shown would he considered as a(t) in the computation of 
the applied landing force F(t) in equation (8). 



Kt) = M a(t) 



(8) 



oscillograph record 
assumed a(t) curve 




Figure 7 



5. EVALUATION OF THE STRESSES 



Having derived the time history of the deformation of 
the structure, the time history of the stresses is obtained 
by a straightforward procedure. In this discussion, stress 
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is used as a general term which refers to shear or moment. 
It is convoniont to consider the total stress as resulting 
from the superposition of the stresses due to the deformation 
in each mode. In this approach, the zero mode which is a 
rigid motion does not contribute any stress. Only the first 
mode, second mode, and so forth, have to he considered. 
From a practical viewpoint it is also important to note that 
the stresses in each mode are conveniently calculated by us- 
ing the inertia forces of the natural oscillation rather 
than hy trying to calculate the strain from the space curva- 
ture of the mo.de shapes. Consider, for example, a pure 
"bending mode as represented in figure (8). 




Figure 8 



The shape of this mode is represented "by a function 
<t>^(x) such that if q x is the tip deflection (Jj^xo^ 
is the deflection at the location x. The "bending moment 
in this mode is 

M(x) - II dS + (l a )(x) <u O) 
d x a 

M is proportional to q x and varies along the span as 
s (l ) 

E I i± — z. LZ-i.. Obviously the process of computing the 

as 3 

second derivative of <J>^ l ^(x) is elaborate and inaccurate. 
It is therefore preferable to derive the bending moment in 
each mode directly from d'Alembert's principle by expressing 
the equilibrium condition which exists in the natural mode 
between the bending moment at station k and the inertia 

2 ( 1 ) 

forces co x m^. <t> y, q. x iue to oach mass m k located out- 
board from x. Similarly the. shear in each mode may be ob- 
tained from the summation of all inertia forces. 
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It is Been that the stress in the jth mode at 

the kth wing station may be expressed as 

= A (J) q (10) 
k k j 

each of these stresses being proportional to the coordinate 
qj of the corresponding mode. The total stress is obtained 

by superposition and is 




(11) 



7. THE DYNAMIC RESPONSE FACTOR Y 



It has been shown that the response of each mode is de- 
fined by , the motion of an equivalent simple oscillator acted 
upon by a generalized force Q(t). The motion of a simple 
oscillator under the action of an arbitrary force F(t) is 
given by an evaluation of Duhamel's integral as shown by 
equation (l). Consider, for oxamplo, the motion of a simple 
oscillator under the influence of the isosceles triangle 
force-time impulse shown in figure 9. 



7 > r r f f — s T~ 



M 



9max 







-<: 


T I * 


t 



Figure 9 



Let q = static displacement caused by Q m __ 

S * 213. X 



T = 2tt/(D = natural period of the oscillator 
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= period of the triangular impulse 

The complete time histories of the motion of the oscillator 
have been evaluated for the isosceles triangular impulse for 
two ratios of Tj/T, and are plotted in figure 10. The re- 
sponse of the oscillator is expressed as a ratio of its ac- 
tual displacement to its static displacement under Q m ax* 




Figure 10 



Tor each period ratio Tj/T there is a maximum value of 
q/q B in the positive direction, and a maximum value in the 
negative direction. These maximum values are designated here 
as dynamic response factors and denoted by "Y with the proper 
sign affixed. There are two values of f associated with 
each period ratio as shown by figure 10. Consequently a curve 
may be drawn showing the variation of dynamic response factor 
with period ratio. Such a curve for the isosceles triangular 
impulse is shown by figure 11. 
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Ij ^ Period of impulse 
Figure 11 Tj Period of oscillator 



The determination of a dynamic response factor-period ratio 
curve may "be accomplished for any arbitrary Q(t) variation 
"by a numerical or analytical evaluation of Duhamel ' s integral. 
However, both processes are quite lengthy, and require plot- 
ting of the time history as shown by figure 10. A mechanical 
analyzer consisting of a torsional pendulum has been devel- 
oped (reference l) which may be used to measure the dynamic 
response factor without recourse to an evaluation of the com- 
plete time history of the motion. By such means, a dynamic 
response factor diagram for any shape of Q(t) curve may be 
evaluated in a relatively short time. 



8. STATISTICAL APPROACH TO THE LANDING PROBLEM 

The methods outlined above are applicable when the time 
history of the externally applied landing impact forces are 
accurately known. Actually, the shape of the landing impact 
force-time curve is apt to vary considerably from one landing 
to the next, and with different operating conditions of the 
airplane. Furthermore, the responses of the various modes 
are sensitive to small variations in the shape of the P(t) 
curve. For these reasons, it is not practical to design for 
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a single mathematical landing, but instead it is desirable 
to employ a statistical approach to the problem. In the 
previous section it was shown that the extreme positions of 
oscillation of a simple oscillator for any one type of ap- 
plied force-time curve may be conveniently represented by 
means of a dynamic response factor-period ratio curve. It 
is possible to consider a large number of shapes of landing 
force-time curves varying from soft to harsh landings, and 
evaluate a dynamic response factor-period ratio curve for 
each of them. These curves may be plotted on the same graph,, 
and an envelope curve may be drawn which bounds all of them. 
This envelope would represent conditions which exceed in 
severity every type of landing which was considered. Such a 
statistical basis may be used to establish design landing 
dynamic response factor envelopes for tho landplano and the 
seaplane. (For example, see fig. 13.) By means of the dy- 
namic response factor envelope, the maximum deflection of 
the structure in each mode during the landing may be evalu- 
ated quite simply. Considering any single mode, it is as- 
sumed that the structure is restrained to deflect in that 
mode only, while the maximum value of the generalized force 
is applied statically. The deflection may be computed under 
these conditions and then multiplied by the dynamic response 
factor in order to obtain the maximum dynamic deflection 
during the landing. 

f v A 6) 

^(max) J ^(static) J M w a J M w a 



It is apparent that the phase relations between the modes 
are lost when an approach of this type is employed. However r 
this is not serious since for design purposes it must be 
assumed that sometime during the life of the airplane the 
phasing between the modes will be such as to produce the 
worst combination of stresses. For this reason, the maximum 
deflections are superposed without regard for phase in order 
to yield the most critical combination. 



9. THE LANDPLANE WING 



In the case of the landplane, the forceB and moments 
are applied to the airplane structure through the shock 
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strut-tire system. These consist of a force applied in an 
axial direction along the shock strut caused by the landing 
reaction component parallel to the strut, and a moment about 
the landing gear attachment point caused by the landing re- 
action component perpendicular to the strut. Two dynamic 
response factors must be determined. One for the axial strut 
force, and one for the moment about the landing gear attach- 
ment point. 

The shook strut axial force. - A study of force-time 
curves for axial strut reactions obtained from drop and 
flight test data shows that they differ considerably with 
pilot technique and landing attitude. A group of six types 
of force-time curves are considered which would represent 
various types of strut characteristics . These are shown in 
figure 12. 




(4) (5) (6) ■ 

Figure 13 

Dynamic response factor curves similar to figure 11 have 
been determined for each of these six curves by meanB of the 
mechanical analyzer (reference l). The pointB which define 
these curves have been plotted on one graph in figure 13, and 
a smooth envelope curve has been drawn which bounds all of 
them. The condition, which the envelope shown in figure 13 
represents, exceeds in severity the condition of the six 
types of landings considered. It is possible that after con- 
sideration of a large number of shapes of axial shock Btrut 
force-time curves taken from drop and flight test records, an 
envelope of this type may well represent conditions which ex- 
ceed in severity every probable landing which would be exper- 
ienced during the normal operation of a land-type airplane. 
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In order to apply the landplane dynamic response factor 
envelope, the impulse period must be known. A plot of 
limited data available at the Bureau of gross weight against 
landing gear vertical impulse period during the first impact 
of landing as dotormined "by flight test is shown on log-log 
paper "by figure 14. This graph has "been determined from ac- 
celerometer records of various types of landings made "by 11 
airplanes of various weights. Bach point represents an av- 
orago of several landings. The equation which fits the 
curve drawn in figure 14 is 

0.14 75 

0.25 I (13) 



{1000} 



*IN 
where 

Tjjj. axial shock strut impulse period in seconds 

W gross weight of airplane in pounds 

The wheel drag force .- The characteristics of the drag 
force on the wheel are not as well known as those of the 
axial strut foroe. The drag force is produced by the spin- 
ning up of the wheel when ground contact is made. If it is 
assumed that the tire is slipping or on the verge of slip- 
ping on the runway at all times during the wheel spin-up 
time, and that the coefficient of sliding friction n is 
constant, the following differential equation may be written. 

D(t) = idT(t) = (14) 

r e 



Integrating once gives 




t)dt (15) 



where 

landing speed, feet per second 

I w moment of inertia of wheel and tire, slug-feet square 

r effective rolling radius of wheel under impact loading, 
feet 



UACA ASH Ho. mO 



19 



2 





3. n 
3 M 


























1 x - 

b 


Li' 


; t 

X 


1 


j 
















I 


A ^ 


1 


< 


X 

p 




c 


i 












1 ° 
1 1 


± 


<£> + 


1 7 




V 
A 


( 


>* 

T 


□ 




































r 


























f — 


























0 




A 
























t 

V 


i 


<? 


I < 

A 


X 

V 


J 














< 


+• 

A 

\ 




: *i 






Refere 
Ml 


nee to 


Figur 


e 12 








x I 


> — i 


j *■ 






(2) A 

(3) + 






7 -l 


o 

> 


-8-] 








(5) 
(6) 


□ o < 






-i-i 


> ( 



















































2 4 6 8 10 12 

Period of irooulqe _ TlN 
Period of jth mode Tj 

Figure 13.- Ifcroamic response factor envelope for landplane axial shock 
strut loads. 



NACA ARB No. 4H10 




(spuooes) Hlj/pof jad esindtnr tboT^A 



NACA ARR No. 4H10 



21 



Tjjj period of drag force impulse, seconds 

If the shock strut axial force is assumed to grow linearly 
with time during the wheel spin-up period, then ' equation 
(15) may "be solved for the maximum value of the drag force 
in terms of the period of the drag force impulse. 



max 



3V l ly 

r e 3T ID 



(16) 



With the assumptions used to write equation (16), the drag 
force-time variation is of the type shown in figure 15. 

D(t) 




Figure 15 

The dynamic response factor curve for the force-time 
relation shown "by figure 15 is given in figure 16. The 
wheel spin-up time is of the order of one-fourth to one- 
fifth of the axial strut impulse period, and the wheel is 
often up to speed before the axial strut force reaches its 
maximum value. 



In the case of modern large aircraft with retractable 
landing gears, the assumption that the landing gear leg is 
rigid in fore and aft bending may be considerably in error. 
Because of this lack of rigidity, there is an additional vi- 
bration mode to "be considered which involves large wheel and 
strut amplitude and very little wing torsional motion. This 
lack of strut rigidity may he particularly troublesome if 
the fundamental fore and aft bending frequency of the land- 
ing gear is coincident with some other mode of the structure, 
and resonance is established. 

In a general consideration of the landplane landing 
problem, drift and one wheel landings which excite antisym- 
metric wing oscillations should he considered as well as sym- 
metric landings. The principles discussed heretofore are 
quite general and apply oqually well in either case. 
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Dynamic response factor curve for 
landplane drag force impulse of 
the type shown l»y figure 15. 




Figure 16. 
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10. THE FLYING BOAT WING 



In the flying boat problem, the forcing impulses are 
transmitted to the wing through the hull structure, and the 
problem may be assumed identical to the landplane problem 
except that the force and moment is applied at the center 
line of the airplane rather than outboard at a landing gear 
station. The problem of the impact of a flying boat hull in 
water is considerably more complex than the landplane impact 
problem. A limited amount of test data is available showing 
time histories of center of gravity vertical and drag ac- 
celerations, and pitching acceleration for flying boats. 
These records indicate that the time history of the landing 
reaction varies widely according to landing attitude, pilot 
technique, condition of seaway, and detailed characteristics 
of hull. Theory on the seaplane impact problem is extensive; 
however, its applicability to the complex hull shapes of 
modern flying boats has not yet been demonstrated, and hence 
little attempt is made to use it for design purposes. Not 
only the force-time relations of the vertical and drag forces 
are necessary, but also their lines of action on the hull 
bottom must be known. Of these necessary items, the charac- 
teristics of the vertical force are more completely known 
than any of the others. Very little general information is 
available concerning the drag force and how the exact line 
of action of the vertical and drag forces vary throughout 
the impact period. Because of this lack of information, it 
is difficult to determine exactly how the impulses are ap- 
plied to the flying boat wing, Furthermore, the influence 
of the moment on the wing vibration is undetermined unless 
the effect of fuselage and tail oscillations are considered. 
Such a consideration is beyond the scope of this discussion. 
A first approximation may be obtained by considering only 
the effect of the vertical force applied at the elastic axis 
of the free wing. In the absence of more complete test data 
it may be assumed that the vertical force on the seaplane 
hull, during the initial landing impact, varies as a half 
cycle of a sine wave. A dynamic response factor envelope 
has been determined for a half cycle of a sine wave impulse 
and is given in figure 17. The variation in the impulse 
period is as wide as the variation in the shape of the im- 
pulse curve, and hence it is not possible -to derive an em- 
pirical relation for the impulse period from test data, as 
was the case with the landplane. In order to apply the dy- 
namic response factor curve, an impulse period must be as- 
sumed, or a value taken from flight test data on an airplane 
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Figure 17.- Dynamic response factor curve for seaplane vertical landing 
force-time curve consisting of & half cycle of a Bine wave. 
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of similar size to the one being investigated. A conserva- 
tive procedure, in any case, would be to choose a period 
ratio so as to yield the most critical combination 

of stresses. 



11. I LLTJ STRAT I VE NUMERIC AX EXAMPLE FOR LAND PLANE WING 



The principles discussed in section 9 » re illustrated 
by a numerical example in which the stresses due to landing 
in a four-engine land type bomber are computed. Mode shapes 
and frequencies of the free wing of this airplane were 
computed and transmitted in a report prepared for the Army 
Air Forces at C.I.T. by M. A. Biot in 19&1. The bending and 
torsional moments are computed at seven wing stations. In 
this example, the wing chord is assumed parallel to the 
ground during the landing, and the landing gear strut is 
assumed to be perpendicular to the wing chord line. The land- 
ing gear strut is assumed to be rigid and rigidly attached to 
the wing. Figure 18 shows the assumed conditions during the 
landing. 



Angular dis- 
placements and 
torques 



Linear displacements 
and forces 



In this particular example, the elastic axis passes di- 
rectly over the landing gear leg as is shown by figure 18. 
The computations may be carried out by table 1 which is self- 
explanatory when reference is made to appendix I. A table 
similar to table 1 is required for each wing mode. The mode 
shapes and frequencies of the free wing are taken from the 
Army Air Forces report prepared at C.I.T. and tabulated in 
tables 2, 3, and U, where the first three wing modes are 
considered. The following additional data are required to 
complete the tables. 
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Gross weight, V 47,200 lb 

Gross weight less weight of landing gear, W L 44,426 l*b 

Moment of inertia of wheel and tire, 1^ 28 slug-ft 2 

Effective rolling radius, r Q 2.08 ft 

Landing speed, V g 124.8 fps 



Prom equation (13), the period of the vertical impact 
force is, 

Ttw = 0.25 f-$-A " =0.25 ' = 0.441 second 

L » \l0O0j 1^1000 j 

The maximum value of the vertical impact force is 

N max = | w L n = I X 44 ' 425n = 22 ,213n 

where n is the ultimate vertical load factor. In this ex- 
ample, it is assumed that n = 4. 

M max = (22, 213) (4) = 88,852 pounds 

It is assumed that the period of the drag force is one-fifth 
that of the vertical force. 

T ID = 5 T IN = 0 *^ 41 = 0.0882 second 
From equation (16), the maximum value of the drag force is 

D m ililiL . (2)(124.8)(28) = 18f290 pQunds 
r Q 3 T ID (2.08) 3 (0.0882) 

From the data computed in columns 5, 6, 13, and 14 in 
tables 2, 3, and 4, wing tending and torsional moments are 
plotted in figures 19 and 20, respectively, for each mode in 
the first and second extreme positions of oscillation. Crit- 
ical values of "bending moment and torsional moment at each 
station are obtained "by adding corresponding ordinates of 
the three bending moment curves on the same side. of the zero 
axis. For example, at the airplane center line, the criti- 
cal negative bending moment is obtained by adding: 
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Inch-pounds 

First mode - positive maximum • -3,940,000 

Second mode - negative maximum ........... -790,000 

Third mode - negative maximum -212, 000 

Maximum negative lending moment at the center line = -4,942,000 

Similarly, for the maximum negative torsional moment at 
the center line, the following are added: 

Inch-pounds 

First mode - positive maximum ............ -1,092,150 

Second mode - positive maximum -73,316 

Third mode - positive maximum -142.525 

Maximum negative torsional moment at the center line = -1 , 307 , 991 

The frequency of the fourth mode is approximately 1350 
epm, and hence its contribution to the stress would he small. 

This is true "because of the tendency for to he reduced 

to small values in the higher modes, and "because of the in- 
hibitive effect which the aerodynamic and structural damping 
has upon the higher modes. It is important to remember that 
the stresses,, shown here must be superposed upon the steady 
stresses produced by the aerodynamic loads on the wing dur- 
ing landing. 

12. ILLUSTRATIVE HUMERI CAL EXAMPLE FOR SEAPLANE WING 



The principles discussed for the seaplane are illus- 
trated by a numerical example. In this example, a four- 
engine land type patrol bomber is considered as being a 
flying boat in order that the same mode shapes and fre- 
quencies may be used. The bending and torsional moments 
are computed at seven wing stations for a vertical load 
factor of 1, and the effect of the drag force is neglected. 
The first three deformation modes of the free wing are 
considered, and the computations are carried out by means 
of table 1. The vertical force on the hull is assumed 
perpendicular to the wing chord during impact. 
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Figure 19.- Iandplane-wing bending moment vs. wing station. 
First three modes for 4G landing. 



iNACA AfiR Ho. 4H10 



33 



O 



o 



100 



Distance from k ■* inches 
300 300 ~" 400 



500 



600 



in 
§ 

i 

•H 

Vt 

O 

w 

o 



•H 

s 



CD 



1.0 



.8 



First mode-negative maximum 



.6 



Second mode-negative maximum 



.2 



Third mode- negative maximum 




Second mode-positive maximum 



Third mode-positive maximum 



— • 8 



-1.0 




-1,2 



First mode-positive maximum 



Figure 20.- Landplane-wing torsional moment vs. wing station. First three 
modes for 40 landing, 
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The following sketch shows the assumed conditions during 
the landing. 




NU) Figure 31 

Por purpose of this example, an impulse period of 0.2 
second, is arbitrarily assumed. Prom the data computed, in 
columns 5, 6, 13, and 14 in tables 5, 6, and 7, wing bending 
and torsional moments are plotted in figures 22 and 23, 
respectively, for each mode in the first and second extreme 
positions of oscillation for a vertical load factor of 4. 
The values plotted in figures 22 and 23 are four times the 
value s shown in columns 5, 6, 13, and 14 in tables 5, 6, 
and 7. 



13. DISCUSSION 



In the usual landing with large aircraft there is a 
period of initial impact followed by a short interval during 
which the airplane is airborne. It is to conditions existing 
during the initial impact and shortly thereafter that the 
present theory of transients is applicable. In the ensuing 
motion during run-out, the wheels are substantially in con- 
tact with the ground or the hull with the water. During this 
latter phase of the motion, the structure has imposed upon it 
a series of alternating loads depending upon the ground con- 
tour or the condition of the seaway. It may be possible that 
a resonance condition during the latter phase of the landing 
may produce stresses more critical than those produced during 
and shortly after the initial impact. 

A preliminary analytical investigation into the effect 
of aerodynamic damping has shown that it may be neglected 
with small error; however, this is subject to experimental 
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Figure 22.- Seaplane-wing bending moment vs. wing station. First 
three modes for 4G landing. 
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Figure 23.- Seaplane-wing torsional moment vs, wing station. First three modes 
for 4G landing. 



NACA ARB. No. 4H10 



40 



check. Damping has the effect of reducing the amplitudes 
and coupling the modes. 

Coupling "between the motion of the structure and the 
external force is neglected in this discussion; however,this 
assumption is also subject to experimental check. 

The methods discussed here are not universally appli- 
cable to all types of transient problems, and should be used 
with discretion. For example, appendix III discusses the 
case of a simple beam subjected to a unit impulse applied at 
the center. In this case, the method is not applicable, 
since reference to equations (III-ll) and (111-13) shows 
that the series describing the moment and shear does not 
converge. However, in practical problems, the applied forces 
do not behave like the impulse type, but rather behave like 
the step type. In this case, reference to equations (III-10) 
and (IH-12) shows that the series describing the moment con- 
verges, whereas the series describing the shear doe6 cot con- 
verge when maximum values are added regardless of sign. In 
the case of the procedure discussed here, convergence of the 
series describing the moment and shear is dependent upon di- 
minishing values of <jj^} This is assured because of the 
characteristic taper of the airplane wing from center line 



side of the wedge in its region of greatest mass per unit 
length. These limitations regarding the convergence of the 
series apply even more so if one attempts to predict the 
local accelerations in the structures. This aspect of the 
problem becomes significant when evaluating the dynamic 
stresses on the attachments of small localized masses. It 
is not possible to formulate at this stage a general rule as 
to how many modes should be taken, and each application of 
the procedure will present a different problem depending 
upon the mode shapes and frequencies of the wing. 

In the landplane case, the position of the landing gear 
leg has an important effect upon the wing dynamic loads dur- 
ing impact. For example, if the leg intersects a nodal line 
for one of the modes, that particular mode is not excited. 
The effect of various landing gear positions on the stresses 
may be readily studied by these methods. 

In studying the present day large aircraft with con- 
ventional wings, it appears as though wing dynamic loads will 



to wing tip, and because 




is measured at the inboard 
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produce torsional moments inboard of the nacelles which, 
aggravated "by the overhanging engines, may "be critical in 
sev-ere impact. The design of the nacelle carry through 
structure in bending may he controlled "by dynamic loads, 
and the possibility of critical wing "bending stresses is 
not precluded in very large aircraft. 

A comparison "between the landplane and flying "boat ex- 
amples givon in sections 12 and 13 shows that the wing dy- 
namic loads in the flying "boat are more severe than in the 
case of the landplane. This is attributed" to two causes: 

(a) In the case of the second mode in the landplane 

example, the drag force actually has a reliev- 
ing effect on the stresses produced by the 
vertical force. 

(b) In the case of all modes (t/jp at the center line 

is greater than <J>^ at station 1 - that is, 

the forces introduced at the landing gear are 
applied nearer to the nodal lines than forces 
introduced through the hull at the airplane 
center line. 



The present work is of a preliminary nature and many 
questions are left for further investigation. It is evident 
that the methods here presented are not restricted to the 
evaluation of landing loads, since it is possible at least 
theoretically to handle in the samo way dynamic loads due to gun 
recoil and "flak." It must be remembered, however, that in 
flight the aerodynamic forces become of primary importance 
and cannot be generally neglected. This is especially true 
in the determination of dynamic loads due to gusts in which 
case the flutter characteristics of the airplane must have a 
preponderant effect. It must also be kept in mind that the 
possibility that the representation of the transient motion 
as a superposition of natural modes is not necessarily the 
best procedure in all cases. Considering the dynamic stresses 
from the standpoint of wave propagation in the elastic system 
might turn out to be a more direct and significant viewpoint 
in the case of high frequency transients. This viewpoint 
also eliminates the convergence difficulties mentioned above 
in connection with the determination of local accelerations. 
Another case where natural modes lose their significance 
is that of resonance between- loosely coupled parts of the 
structure, in which case the vibrational energy at one loca- 
tion is gradually transferred to another. 
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Examples shown in the present work are limited to the 
wing structure under the assumptions of a symmetric landing. 
Such a landing condition is exceptional. For a landplane 
the degree to which the unsymmetric modes are excited by an 
unsymmetric landing depends a great deal on the time inter- 
val between the instant at which the left and right wheel 
enter in contaet with the ground. Statistical data on this 
time interval can only "be obtained by flight testing. No 
example has been presented of an application of the proced- 
ure to the evaluation of landing loads in the fuselage and 
tail. However, the same methods are directly applicable to 
this case, provided the natural modes of the fuselage and 
tail have been determined. Data obtained during landing 
tests of flying boat 6 have shown that modes of the fuselage 
and tail are excited and result in a "tail whip" effect 
causing considerable dynamic overstress in the tail and aft 
portion of fuselage. 



Bureau of Aeronautics, 

TJ . S. Navy Department, 

Washington, D , C., August 10, 1944. 
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APPENDIX I 



GENERAL MATHEMATICAL THEORY OF TRANSIENTS IN AN 
UNDAMPED ELASTIC STRUCTURE 



The general transient theory of linear systems with 
lumped or distributed parameters is well known and has "been 
extensively developed in the case of electrical network 
theory (references 2 and 3). The problem of transients in 
airplane structures is identical in its mathematical form. 
In an elastic system with distributed parameters there is a 
space as well as a time variation in the variables. The 
problem may be considered with two viewpoints. The motion 
may be considered to be made up of a series of traveling 
waves, or it may be considered to be made up of a superposi- 
tion of natural oscillations, in which case to be rigorous, 
an infinite number is required. The airplane structures 
problem is treated here from the standpoint of a superposi- 
tion of natural oscillations. 

In a transient problem of this type where maximum val- 
ues occur very soon after the motion starts, the effect of 
damping may be justifiably neglected. The motion of an un- 
damped elastic system may be shown to be composed of a super- 
position of normal modes which are orthogonal. The airplane 
structure vibrates in a series of normal modes when excited 
by a random impulse as is the case of any elastic system. 
These normal modes are each characterized by a certain mode 
shape and a certain frequency. For the airplane they are 
composed of coupled oscillations of the wing, fuselage, and 
empennage system. These calculated mode shapes and fre- 
quencies may be obtained from a ground vibration survey of 
the airplane. 

If each normal mode shape is considered to be repre- 
sented by the space function * ' the displacement of any 
point of the structure at any time may be written as 



n 




(1-1) 



i = o 
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where n modes are- considered, and where the q^ terms are 
regarded as generalised coordinates. 

The kinetic energy of the structure is 

n 



<fr(i) qi j 3 dm (1-2) 

o 



Considering the jth generalized coordinate, and the orthog- 
anality condition, the following equation may be written 



j_ /dT 
dt 



where 

I a 



"i"/[« (J, I 



dm 



The potential energy of the structure is 

n 

U = I V *i ">i* (1-4) 
i = o 

Considering the jth generalized coordinate, write 

£H- = Mj Wj » q, (1-5) 

By Lagrange's equation 

M d 4j + M j ^ j 3 4j = *.J (I ~ 6) 



where is the generalized force corresponding to the jth 

mode, and is evaluated "by virtual work principles. 
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If 6Wj is the virtual work produced when all the ex- 
ternal forces are allowed to move through displacements cor- 
responding to a virtual displacement 8qj, then the gener- 
alized force Q,j 1b defined "by 

Ij-TT* (1-7) 

The complete motion of the structure is then defined by a 
series of differential equations of the form 



where the form of the right-hand side is dependent upon the 
character of the applied forces. 

In general, the aerodynamic applied forces on an airplane 
structure vary with deflection, velocity, and acceleration of 
the structure, and the landing reactions vary with time in a 
manner which is determined hy experiment. In the case where 
the external forces are landing reactions assumed to be givon 
functions of time, the equations governing the response are 

«j + »/ *j - ^Sp (I - 8) 

This is the differential equation for the undamped mo- 
tion of a simple oscillator of mass Mj and natural frequency 
U)j, which is under the influence of an arbitrary forcing 
impulse Q j(t ) , 

If Q^(t) is a unit step function l(t), the response 

q. which is called the indicial admittance (references 3 
J 

and 4) is 



A(t) = i (1 - cos u>4 t) l(t) (1-9) 
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The response q.j to any arbitrary forcing impulse Qj( t) 
may be written by the superposition theorem as 

t 

q . a — - — f Q,At) sin CO, (t - t) dT (1-10) 

T is a variable of integration. 

When this integration with respect to t is carried out 
between the limits 0 and t, a function of time results 
which is the time history of the deformation of the jth mode 
The stress at any point in the structure in the jth mode is 
proportional to the deformation of the jth mode, 

(3) 1 

s^ 5) = q. = — f Q.(t) sin CO. (t - T)dT (l-ll) 

When th© constant A^ is properly chosen, equation (I-ll) 
yields the stress time history of some particular point in 
the structure caused by the deformation of the jth mode. 



.(J) 
' k 

may be written as 



The stress s^ ; in the jth mode at the kth wing station 

is 



The total stress at the kth wing etation for n modes is 
obtained by superposition as 



s k = £ A-' q, (1-13) 



Note on the Computation of Mj for the Wing 
From equation (1-3) 



dm 
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For the case of the wing, the mode shape is conveniently de- 
scribed by considering the wing deformation to be made up of 
a bending of the elastic axis, and a twisting about the 
elastic axis. Considering the wing to be divided into k 
spanwise stations, the normal function describing any point 
on the chord of the kth station is (see fig. I-l) 



+ x a 



(J) 



k 



(1-14) 




Figure I-l 

In figure I-l positive bending deflections are downward and 
positive pitching deflections are stalling. 

Using equations (1-3) and (1-12), the following equa- 
tions may be written 

m 3 = y, {[^ j) ] s ^ * [ a * j> r ** + ^ (i - i5) 



Note on the Computation of the Generalized Force Qj 
From equation (l-7) 



W 
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If the external force F(t) is a landing reaction consid- 
ered to be applied at one point to the wing structure, the 
virtual work may he written as 

SW * T(t) 6 zy 

where zy is the value of z in the direction of the ap- 
plied force evaluated at the point F, the point of appli- 
cation of the landing force . 

From equation (I-l) 

. q x + <t> (3) *,+ ... * (n) q n 



hence 



r <o j*) a*) -\ 

8W = F(t) <<> p 8q x + $ P 6q s + . . . * f 6q n | 



where 0^*1 is the normal function evaluated at the point 
F 

of application of the force F. 
Then 

Q J3L „ F(t) (1-16) 

The factor is a measure of the contribution of 

F 

the external force to the generalized force in the Jth mode. 
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APPENDIX II 

ANALYSIS OF THE MOTION OF A SYSTEM CONSISTING OF TWO SPRING- 
CONNECTED MASSES EXCITED BY AN ARBITRARY FORCE 



fx 2 












k 
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7 — r ~r~ 


m 
L r- r 




/ f ' r / t 
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7 / 





F(t) 



Figure II-X 



It is shown "below that the motions of the mas6 m are 
given hy a superposition of the motion of the center of 
gravity of the system, and the motion of a simple oscillator 
which will be defined. Denote by xj., X 2 the coordinates 
of the masses M and m. 

The differential equat i ons d e fining the motion are: 

Mx x + k( Xl - x 2 ) = F(t) (II-1) 
mx 3 - k(x 1 - x 2 ) = 0 (II~2) 

The motion of the center of gravity of the system is de- 
termined first. Adding equations (il-l) and (II-2) gives 

M Xl + mx s = F(t ) 

Let q 0 be the displacement of the center of gravity of 
the system 

rnx 2 + M Xl = (M + m) q 0 

and 

mx 2 + Mx\ - ( M + m)q" 
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Hence 

(M + m) q Q = y(t) 

is the differential equation defining the motion of the cen- 
ter of gravity of the system. 

The motion of the mass m relative to the center of 
gravity of the system is determined next. 

Multiplying equations (II-l) and (II-2) by m and M, 
respectively, and subtracting, the following equation is ob- 
tained 

Mm(x l - x s ) + k(M + m)(x l - x a ) = ml(t) 
and substituting 

4 X = rr*— <*i - x 2 ) 
M + m 

which represents the motion of m relative to the center of 
gravity; the differential equation may be rewritten as 

■ (»♦«) 5l+ . k (aii)\.*f(t) (ii-4) 

Let 

H 1S { (m + JO 

M 

- jf »(♦) 

and write 

Mi qi + Mi^ 3 q x = t^CO (II-5) 

This is the differential equation defining the motion of m 
relative to the center of gravity of the system. It is also 
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the differential equation for a simple oscillator. (See fig. 
II-2.) 







— *1 


Ml 








— / ■ r / -r ■/ / / 




* r r "f 


Figure II-2 





Qx(t) 



The actual displacements * x and x 3 may te derived by 
solving for x 1 , X 2 . 

From the relations 

mx 3 + = (M + m) q 0 

may be found 

m 

x i = Qo + ^ 1i 

The motion may he considered as the superposition of two 
configurations, one defined by q Q , a rigid motion of the 

system, the other defined by q x represents a configuration 
in which the center of gravity remains fixed while the masses 
M and m move in opposite phase with amplitudes inversely 
proportional to their masses. The equivalent system of fig- 
ure 11-2 represents the motion in the latter configuration. 
It will be noted that the generalized mass Mi may be derived 
quite simply by considering th-e kinetic energy T in the 
corresponding configuration 
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or 



T = I J»(m + M) q ' 
2 M Hl 



Equating this to the kinetic energy of the equivalent oscil- 
lation 

T = I M x q\ 2 
2 



there is found 



M, = £ (m + M) 
M 



Similarly Q x (t) may be derived by equating the work done by 
F(t) in the corresponding configuration and the work done by 
in the equivalent system. 



Q x (t) <u - F(t) £ q a 



Hence 



/ 



<4x(t) = g l(t) 



APPENDIX III 



TRANSIENTS IN A PRISMATIC BEAM SIMPLY SUPPORTED AT THE ENDS 

with a force applied at the center 

,F(t) 



1 



Figure III-l 
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The natural mode shapes and frequencies of a simply sup- 
ported team with constant cross section may be shown to "be: 

= sin 125 

v 



0), = iV AHl (III-2) 
VoAl 4 



wher e 

B modulus of elasticity 
I moment of inertia 
p mass per unit volume 
A cr o s s- sect i onal area 
I length 

Prom equation (i-l) of appendix I, the displacement of any 4 
point on the beam may be written as 

00 CO 

y = £ 0<O qi = £ 8in iSi qi (III-3) 
i=o i=i 



From equation (1-10) of appendix I, the generalized coordinate 
q^ is expressed by 

qi = 77-77- / M T > Bind), (t - T) dT (III-4) 
i i 

where j 

Mj _ = J [o(i)] S dm = pA f sin 2 ±2* dx - Ig± (III-6) 
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and 

Q t (T) « P(T) = sin AH P(T) (III-6) 

Substituting equations (III-4), (III~5), and (III-6) into 
equation (III-3) yields, 



Z, t 
ill sin T Bijl l fl P(T) 8in B i (t ~ ^a^ 111 - 7 ) 

K .1=1 o 



When P(t) is a unit step, the response is: 



7 = -2- 
IpA 



■3D 

) _i_ sin iH sin l££ (1 - cos (» ± t) (III-8) 
£-j W x s 2 1 1 



When F(t) is a unit impulse the response is obtained 

by a time differentiation of equation (III-8). The response 
to a unit impulse ie: 



CO 

2 



) JL B i n iEL sin I2i sin <»! t (III-9) 



7 " IpA 

i=i 

Since M = -EI the moment in the beam at any time after 

being subjected to a unit step force is, 

CO 

M a 21 y 1 Bin JLE Bin in* (l - C08 (i), t) (111-10) 
l = i 

The moment in the beam at any time after being subjected to a 
unit impulse is, 

( 1 ^A unit impulse is an infinite force applied during a 
vanishingly small interval of time in such a way that the 
time integral / ^(T) at = 1 (reference 4). 
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00 



M = | Vsin If sin ^ sin w 4 t (111-11) 

i = i 

The shear in the "beam at any time after being subjected to a 
unit step force is, 



M = 1 V 1 6i n^co 8 lH(i- cos U), t) (111-12) 
i = i 

The shear in the "beam at any time after "being subjected to a 
unit impulse is, 

00 

dM 2tt /EI V 1 . . in Ittx . . / TTT ,„\ 

— = S — / i sin — cob sin to* t (111-13) 

3x W pi L 3 , I 



A study of equations (III~10), (ill^ll), (111-12), and (111-13) 
indicates the shortcomings of the procedure when adding con- 
tributions of each of the modes to moment and shear for a 
simple prismatic beam with a force at the center. This is 
illustrated in table III-l. Limitati ons of a similar nature 
are encountered in all problems where the motion is described 
as a superposition of modes. 



TABLE III-l 



UNIT STEP 


UNIT IMPULSE 




Percent of 1st mode 




Percent of 1st mode 


Moment 




Moment 




1st mode 


100 


1st mode 


100 


2nd mode 


11.1 


2nd mode 


100 


3rd mode 


4 


3rd mode 


100 


4th mode 


2.04 


4th mode 


100 


Shear 




Shear 




1st mode 


100 


1st mode 


100 


2nd mode 


33 .3 


2nd mode 


300 


3rd mode 


20 


3rd mode 


500 


4th mode 


14 .3 


4th mode 


700 
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